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We provide an intuitive interpretation of the optical Stern-Gerlach effect (OSGE) in the dressed- 
state point of view. We also analyze the effect of atomic damping in an experiment on the OSGE. 
We show that the atomic damping also causes the wave packet splitting, in a non-mechanical fashion, 
as opposed to the coherent process that is mechanical. 
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I. INTRODUCTION 



The mechanical interaction of light and matter has 
been comprehensively studied since the pioneering work 
of Kapitza and Dirac in 1933 [0. Much work has been 
done on the atomic beam deflection, diffraction, refrac- 
tion or interference H-j^ by a standing-wave field, classi- 
cal or quantized, and recent achievement of optical cool- 
ing and trapping of neutral atoms |Q using the light pres- 
sure force may be said to have a root in these themes. 
One of the interesting issues that came up in the theme 
of atom-field mechanical interaction is the so-called op- 
tical Stern-Gerlach effect (OSGE) It has been 
explored at various levels of sophistication since mid sev- 
enties when it was suggested that the trajectory of a two- 
state atom interacting with an optical field gradient can 
be split into two paths, each path containing atoms in 
one of the two orthogonal dressed states, under certain 
circumstances It is thus an optical analogy of the 
well-known magnetic Stern-Gerlach effect (MSGE) 
in which the trajectory of a spin-i particle is split into 
two paths in a magnetic field gradient. 

The OSGE was experimentally demonstrated in 1992 
by Sleator and others in the near-infrared (A ~ 1 ^m) 
with metastable helium atoms (He*) [Q. The experi- 
mental result seems to be in agreement with the semi- 
classical theories [^-10|, but later it was also pointed out 
that the OSGE may show some additional peculiarities, 
sensitively depending on the quantum nature of the field 
I p^ . Not only incorporating the quantum nature of light 
and the atomic center-of-mass motion in the picture, re- 
cent works also take into consideration the finite spatial 
extent of the atomic wave packet rather than 

treating an atom as a point mass |^ or a plane wave 
in space. Particularly, the first reference in Jlif presents 
quite an elegant operator method in such a full quantum- 
physical description of the OSGE. 

Noting, however, that the theories originally envisaged 
the OSGE as the coherent processes dealing with ideal 
lossless systems, mostly in the linear regime, we explored 
the role of the decoherent processes in the OSGE, par- 
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ticularly in the nonlinear regime, which may be signifi- 
cant in the optical frequencies. For the systems free of 
damping, decent analytic treatments may be available, 
but when the system is open to its environment so that 
the coherent dynamical evolution of the system is fre- 
quently interrupted by the discontinuous processes such 
as atomic spontaneous decays or cavity emissions, etc., 
analytic approaches may not be always possible. There 
are a number of theoretical methods in dealing with such 
open quantum systems, but we will resort to the formal- 
ism of the quantum trajectory theory (QTT) [ pT| partic- 
ularly because the theory provides most direct and in- 
tuitive interpretations on the behavior of the system we 
look into. In particular, we will show in the viewpoint of 
QTT that the atomic damping can also cause the wave 
packet splitting, and that in a non-mechanical fashion, 
as opposed to the coherent process that is mechanical. 

In Sect.||, we first provide our simple and intuitive il- 
lustration on the OSGE, based on the dressed-statc pic- 
ture and in Sect. [II, we discuss the effect of damping 
in the language of QTT. In Sect. IV, we perform a numer- 
ical simulation of a typical experiment and compare the 
outcomes of such a realistic experiment with the predic- 
tions made for the ideal systems free of damping. Sect. ^ 
summarizes this work. 



II. THE OPTICAL STERN-GERLACH EFFECT 



A. The wave packet "pulsation" 

Consider a simple model in which a Gaussian wave 
packet of a two-state atom initially in its ground state 
is placed on a node of a standing wave cavity, on the 
cavity axis, as illustrated in FIG.nl In order to quickly 
appreciate the essential physics, let us assume that the 
cavity field is initially in a single-quantum Fock state, and 
the atom and the cavity are on resonance. We will neglect 
damping for the time being. The interaction Hamiltonian 
is then simply written as 

H ^ — +ifig{x) {a-a^ - a+a) (1) 

where p is the atomic momentum operator conjugate to 
the position x on the cavity axis; M, the atomic mass; 
and g{x), the atom-field coupling strength at position x 
such that 

g{x) = gosinkx, (2) 

go being a constant and k the field wave number. at(a) 
is the creation(annihilation) operator for the field, and 
a+ = \e){g\ — (c-)^ is the atomic excitation operator, 
with e{g) labehng the atomic excited(ground) state. 

We will focus only on a spatially well- localized — i.e., 
a particle-like — wave packet such that its initial width 
Aa; is much smaller than the wavelength A of the cavity 
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field. We choose Ax = A/407r which corresponds to the 
width in the momentum space Ap = lOhk. Due to the 
finite spatial extent of the wave packet, the atom-field 
coupling strength at the node is nonzero, although very 
small compared to its peak value go, and the interaction 
slowly commences. Note that the initial effective cou- 
pling strength is given by an overlap integral of the mode 
function and the atomic spatial probability distribution. 



geS = / dxt/j^x) \g{x)\ 
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~ 2V2^ (^^^ go ^ 0.0450 (3) 

where "00 i^) is the initial atomic wave packet and we put 
s'mkx ~ kx as Ax <C A. Thus the effective coupling 
strength would be only about 4 % of 50 initially. 

Let us first present our numerical results on the evo- 
lution of the atomic wave packet in FIGs.||(a) and (b). 
Note that since the atomic mass M, field wave number k 
and atom-field coupling constant go appear in the Hamil- 
tonian at the same time, one needs to specify these quan- 
tities somehow, but the only system-dependent parame- 
ter is the dimensionless factor /i — Tik^ /2Mgo, which 
appears from the Schrodinger equation when we rescale 
time in the unit of go^ and the momentum in Hk. [Note 
that = k{hk/M)/2go = kvi/2go is the ratio of the 
Doppler shift by the single photon recoil to the single 
photon Rabi frequency.] Let us take fi = 1.7 x lO""* 
(with M = 100 a.u., A = 600 nm, and 50 = 200 MHz) 
for a typical atom interacting with an optical field in the 
strong coupling regime. 

The figures show that the initial Gaussian packet 
evolves into a double-peaked distribution and then back 
to its original shape, as time proceeds. The inset is a 
ground-level view of the packet in the coordinate space 
along with the cos(A:a;) curve, which clearly shows that 
the two split bumps make turns exactly at the neighbor- 
ing nodes. Let us first note that the two figures are me- 
chanically consistent in the following sense. Let t = got. 
In the momentum space, the position p of one of the 
split peaks grows roughly linearly in time in the fashion, 
p{t) ~ O.Qhkr until r ~ 100, and then bounces back to 
the original position. So, the position i of a peak after 
one cycle of motion in the momentum space is roughly 

plOa/go - 

X ~ 2 y j^dt - (2.9 X 10V)A = 0.49A - 0.5A (4) 

which is in a good agreement with FIG.^(b). Although 
this wave packet motion can be said to have essentially 
the quantum-mechanical origin in the sense that x van- 
ishes if h does, the motion conforms at the same time 
to the classical-mechanics. We will get back to this ap- 
parently trivial observation when we include damping in 
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the system in Sect, y, and will point out some interesting 
features in connection to this issue. 

Thus we have seen an exotic phenomenon of wave 
packet pulsation on a node. Note that the analytic ap- 
proaches such as JT4[ based upon the linear approxima- 
tion did not predict this pulsation. The descriptions are 
mostly limited to the wave packet splitting in an early 
stage of the atom-field interaction. If the atoms are let 
to leave the interaction region at some point of time, the 
emerging wave packet will have a pair of distinguishably 
separated probability density peaks. This phenomenon 
of wave packet splitting is the very optical Stern-Gerlach 
effect mentioned in the introduction. In the following, let 
us provide a simple intuitive account for the physical ori- 
gin of the wave packet pulsation that we have just seen, 
in the viewpoint of the dressed state picture IQ . 

B. Dressed state picture, on-resonance 

Dressed states are the eigenstates of the atom-cavity 
combined system which form a manifold of infinite hi- 
erarchy called the Jaynes-Cummings ladder. Since only 
one quantum is in the system, we need only to consider 
up to the first couplet such that 

|±>^(l/V2)(|0,e)±z|l,g)) (5) 

which have the Rabi-split eigenenergies given by 

E±{x) = h[ujo ± go sin(fca;)] (6) 

where and 1 count the number of quanta in the field. 
FIG.|^(a) shows the spatial variation of E±{x) on the 
cavity axis around a node. In the figure, one can imme- 
diately see that if the system is in state |-|-), it will start 
to slide down the slope of its energy curve and move to 
the left, and in |— ), to the right. However, as one can 
easily show, the state function of the system at an arbi- 
trary time turns out to be always a 50-50 superposition 
of |-|-) and |— ) when the system starts out in an energy 
eigenstate, e.g., Thus, half of the wave packet sit- 

ting on a node is pushed to the left and the other half, 
to the right, resulting in the splitting of the packet into 
two pieces. When the two split bumps of wave packets 
hit the turning points at the neighboring nodes, they re- 
turn just like a pair of classical particles. This is the 
essential explanation as to why we have such pulsation 
of wave packet as in FIGs.^(a) and (b). We can do a 
very simple Newtonian mechanics here again. One of 
the two components of the wave packet will receive from 
the energy slope a force of size F = hkgo cos kx. In the 
early stage of time where linear approximation is valid, 
F ^ hkgo and the peak will gain the momentum grow- 
ing as p ~ hkr initially. But the overall behavior is 
rather close to our eye-estimation of P{t) ^ 0.9?ifcr in 
Eq. (U). One can easily show that the period tq of the 
wave packet pulsation in the coordinate space is exactly 
given by TO = « 400. Cf. FIG.|(b). 
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In this fashion, the dressed-energy gradient created by 
the nonuniform spatial structure of the cavity field plays 
as a dressed-state selector, splitting the atomic wave 
packet into two distinct groups of orthogonal dressed 
states, 1+) and |— ), which is the very OSGE. Let us 
add that, on an antinode, on the other hand, the |+) 
component of the wave packet tends to be pulled out 
to both sides of the antinode while the |— ) component 
tends to be squashed toward the antinode as illustrated 
in FIG. ||(b), giving the net effect of distorting the initial 
Gaussian packet into a wing-broadened, peak-sharpened 
feature. 



C. Off-resonant case 

In the actual experiments, however, it is not practi- 
cally possible to keep the exact resonance indefinitely. 
So let us consider the case in which the atom and cav- 
ity are slightly detuned. Let S = ujc — t^A denotes the 
detuning, where ix'yi(c) is the atomic transition(cavity 
resonance) frequency. The atom-cavity detuning causes 
anti-crossing of the energy levels which makes the picture 
slightly complicated. The first couplet is now designated 

by 

cos6'|0,e) +i sin 6*11, g) (7a) 
= sin6l|0,e) -icos6l|l,g) (7b) 

where 

tan e ^ [5/2g{x)] + ^[6/2g{x)f + 1 (8) 

with the associated eigenenergies given by 

Ef {x) ^fi[u^ + {512)] ± h^g^{x) + {5/2f (9) 

as depicted in FIG.|4| One has to be careful about 
association of the eigenenergies to the eigenstates. In the 
region where X = 5/2g{x) is positive, \d±) are associ- 
ated with E^, but in the negative X region, these are 
associated with E^. This is simply because sin 6*1 ±x = 
cos9\^x (due to the mutually-symmetric behavior of 
sin 6' and cos 6* around 9 = 7r/4), and the eigenenergies 
are determined by the mod-squares of the coefficients, 
i.e., 

E+ = |cos0|^ huja + |sin6'|^ Ticuc for X > 

= |sin6i|^?ia;a + |cos6'|^?iWc for X < 0, (10) 

etc. It is sufficient to consider the solution of Eq. (||) in 
the range < 29 < n. Let us first consider two limiting 
cases. Firstly, in the regions relatively far from the node 
where |X| <C 1, we have 9 — > 7r/4, and therefore \d±) 
respectively approximate to |±). Thus in these regions, 
the wave packet behaves just like in the case of exact- 
resonance discussed above, as expected. Note that, in 
the left hand side of the node, |c?+) has the lower energy 
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E^(x) which approaches E+{x) in Eq. (^; and the 
higher energy (x) « E- (x) , whereas in the right hand 
side, the levels are reversed. Secondly, in the region near 
the node on the other hand, we have d/2g(x) — *■ ±00. 
When g{x) approaches the node from the positive side, 
9 — > 7r/2, and when it approaches from the negative side, 
6* — > 0. Thus the behavior of the first couplet is summa- 
rized as 



g{x) < g{x) > 



Et: M->^-|f,.g) 
E^ : \d+)^\0,e) 



M+>^|1,5) 



(11) 



What it means is the following: if the system is in \l,g), 
for instance, the system is placed in a potential well de- 
fined by E^{x) in this region, and when it is in |0, e), 
it is sitting on a potential hill defined by E^(x). Thus 
if the system starts from the portion of the wave 

packet in this region tends to be trapped in the poten- 
tial well of E^ {x) while the wings of the packet outside 
this central region will slide down the slopes of E^{x) 
on both sides of the node, since outside the region, \l,g) 
is a combination of \d±) « |±) as aforementioned. Note 
that the overall state of the system at arbitrary time 
is no longer a 50-50 combination of \d±). Nevertheless, 
the packet motion occurs always in a symmetric fashion 
about the node, as a matter of course. To see this, write 
the internal state as — Ci|0,e) -I- C2\l,g) with com- 
plex coefficients Ci and C2. Then, at a position x = xi 
(with respect to the origin at a node), it is written as 
li/j) = A+\d+) + where A+ = Ci cos9 - iC2 sin 6* 

and A_ — Cisin^ -I- 16*2 cos0, whereas at a; = ~xi, 
it is given by lip) — B+\d+) + B-\d-) where — 
Ci s\Ti9 — iC2 cos 9 and i?_ = Ci cos 9 + iC2 sin0. But 
since \A±\ — the portions of the wave packet in the 

upper energy state on both sides are the same. Thus the 
symmetry arises around the node. The ratio of the por- 
tions in the upper and the lower energy states, however, 
are changing in time. As a result, we will have the initial 
single-bumped wave packet evolving into a irip/e-peaked 
structure on a node in the presence of finite atom-cavity 
detuning instead of the double-peaked shape in case of 
the exact resonance. 

Our expectation is nicely confirmed by a numerical cal- 
culation as shown in FIG. ^ where we take S/g^ = 0.1, for 
instance, with other parameters the same as before. Of 
course, due to the small effective coupling constant, the 
system slowly evolves from |1, (?) to |0, e), and the portion 
of the wave packet in the central region will also grad- 
ually migrates into both sides of the node. The volume 
of the central region is confined by \{2gi^/5) sin(fca;)| <C 1 
or \x\/X < (l/47r)((5/go)- For Ax/A = l/407r, we have 
|a;|/Ax < 10((5/5o)- If 5/ go ~ 0.01, then |x|/Aa; < 0.1, 
which implies, though loosely, that by far the most part 
of the wave packet resides outside this region and the 
central peak will not be significant if S/g^ < 0.01. 

Let us leave a brief note regarding the dependence on 
the field state. If the field state includes a nonzero ampli- 
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tude of the vacuum state |0) as in the case of a coherent 
state, this portion of field will tend to leave the atomic 
wave packet intact in time. So, there will be always some 
portion of the packet standing still at its original position. 
This will show up as a peak at the center between the two 
separate peaks on both sides, making the overall look of 
the wave packet somewhat similar to the case of nonzero 
detuning aforementioned. This fact is most easily seen in 
the Schrodinger picture. We will get back to this point 
at an appropriate place in the next section since we will 
explicitly deal with the Schrodinger equations there. 



III. THE EFFECT OF DAMPING 



The picture given so far in the dressed-state formalism 
provided a qualitative but intuitive understanding of the 
phenomenon. Nevertheless, the real systems are after all 
open quantum systems subject to damping, and the dy- 
namical effects of such decoherent processes are yet to be 
discussed. The open systems can be dealt with by a num- 
ber of ways, e.g., the master equations the quan- 
tum Langevin equations ]l9| , or the Fokkcr-Planck equa- 
tions 1 19 2^ in appropriate situations. However, since 



we are dealing with the dynamics of a single atomic wave 
packet instead of a statistical ensemble, we find it most 
appropriate to look into the system in the viewpoint of 
the quantum trajectory theory (QTT) |l^. The theory 
provides most direct and intuitive descriptions as to the 
dynamics of such open quantum systems. Furthermore, 
as we will show below, the QTT suggests an interest- 
ing phenomenon that the atomic spontaneous emissions 
in fact catastrophically split the wave packet into two 
or more peaks when the atom is located at or moving 
through a node of the field. 

A quantum trajectory consists of the piece- wise con- 
tinuous coherent evolution and the discontinuous quan- 
tum jumps, such as the atomic and cavity decays, which 
randomly interrupt the coherent evolution. Let Ti (i = 
1,2,...) denote the various decay rates in the system, and 
Cj the corresponding "collapse" operators, e.g., cf- for 
the atomic emission and a for the cavity transmission. 
Then the coherent evolution of the system is given by 
the following Schrodinger-like equation 



Thus the Schrodinger process is slightly modified by the 
damping terms. Because of the damping terms, it is not 
a norm-conserving process, and \^) represents an unnor- 
malized wave function. When the collapse probability 

is greater than a random number R G [0, 1) taken during 
the time interval [t,t -\- dt), a quantum jump occurs in 
the fashion 
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(14) 



Otherwise, the system dynamics follow E g . (|l^ . For fur- 
ther details of the theory, see, e.g., Ref. 17i. 

Now let us return to the original problem in which 
an atomic wave packet in its ground state is placed at 
a node of the cavity field mode, on the axis of the cav- 
ity. The cavity can contain any photonic state of light, 
and let us neglect the cavity damping (k = 0) for the 
time being in order to focus on the effect of the atomic 
damping, and assume a very massive atom so that the 
kinetic energy term can be neglected (Raman-Nath ap- 
proximation). For algebraic convenience, let us write 
g{x) = gQCOs{kx), and assume a very massive atom so 
that the kinetic energy term can be neglected (Raman- 
Nath approximation). Then the Hamiltonian in the in- 
teraction picture becomes as simple as 



H = ihgocos{kx) (cr_a^ — cr+a) . 
We expand the wave function in the fashion. 



(15) 



I*) 



where n is the field quantum number and q the atomic 
momentum quantum number scaled in unit of hk. The 
coherent evolution is then given by the following set of 
dynamical equations: 



-Cg{n,q)\n,q,g)] + Cg(0, g)|0, 



(16) 



-^Ce{n,q) 



[Cg{n,q-l)^Cg{n,q + l)\ - ^Ce{n,q) {n > 1) 



-Cgin, q) = +^ [Ce{n, g - 1) + Ce{n, q + 1)] 



{n > 0) 



(17a) 
(17b) 



where we used the relations cos{kx) = {l/2)[ex-p{ikx) + 
exp(— jfca;)] and exp{±ikx)\q) = |g ± 1) with r — got, 
and 7 — jA/go, the atomic decay rate scaled in g^. So 
the n-th Jaynes-Cummings couplet is decoupled in dy- 
namics from the rest of the infinite hierarchy of couplets. 
Furthermore, it is trivially seen that the vacuum state 
component Cg(0, q) does not change in time. 

The state of the cavity field can be always ex- 
pressed as a linear combination of number states, i.e.. 



(18) 



Let the initial atomic position x — ^qX and the momen- 
tum spread Ap — {Aq)hk. Then the initial wave packet 
is given by 



/oo 
dg^exp 
^ — n 



4(Ag)2 



exp {-2TTiq^o) F{n)\n, q, g) 



(19) 



where F{n) denotes the normalized coefficients such that 

1 



F{n) 



fin) 



(20) 
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The initial values are therefore 



C,"Kg)=0 (n>l) 



C"Jn,q)^F{n)exp 



(21a) 

exp {-2Triq^o) . (21b) 



_4(Ag)2_ 

Now at a node, = 1/4, and Eq. ( ^lh| ) is written as 

C°{n, q) = F(n)e-'3«'e-*^« (22) 

where /3 = l/4(Ag)^, being the width-parameter of the 
initia l wa ve packet. With this initial value, we solve 
Eqs. ( |l7a| ) and ( |l7b| ) in the Eulcrian fashion such that 

C:^;^ (n, q) = Cil-'^^^in, q) - drC^^^'^'^^n, q) (23) 

where time r is divided by small discrete steps dr so that 
r = vdr with an integer v. We obtain the behavior of 
the system at an early time r given by 



Cl{n,q)^-zC°Jn,q) 



Cl{n,q)^C"(n,q)\l + ? ^^"^ 



,4/3(g-l) _ 2 + e~^''('?+^) 



(n>l) 

(71 > 0) 



(24a) 
(24b) 



keeping terms up to the order of (dr)^. Note that by 
comparing with the numerical solutions, we find these 
expressions are valid up to r ~ 5. [At this moment, let 
us take a brief pause and look at the dependence of the 
shape of the split wave packet on the field state men- 
tioned in the last section. The momentum distribution 
V{q) at an arbitrary time is given by 



r{q)^Y.[\Cl{n,q)\' + \Clin,q)f 



K{0,q)\ 



l/(0)P 



: exp 



2(Ag)2 



(25) 



So it is clearly seen that when the field state is contami- 
nated by a nonzero vacuum state component, a miniature 
of the original Gaussian packet will always be seen be- 
tween the separated wings.] If the damping is weak, there 
will be little difference in the coherent evolution in the 
early times during which exp(— 7r/2) « 1. But what it 
does in the jump process is remarkable. Since now there 
is a small probability built up that the atom is in the 
excited state, there is a probability that the atom can 
decay. Let us just assume that such a decay happened 
now. The wave packet then undergoes a quantum jump 
such that 



l^f) — > exp(-ifca;77)cr_|*) (26) 

where exp(— jfca;77) describes the momentum recoil that 
the atom gets from the emission, projected on the a;- 
axis, and is a random number in the range [—1,1]. 
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1 4") represents an unnormalized wave function. In the 
coefScients, this process is written as 



Cl{n,q) 







Cl{n + l,q + 7^) oc F{n+l) 



Since then, the entire wave packet is determined by 
Ggln, q), yielding the probability distribution in the mo- 
mentum space such that 



It is simply the square of the difference of two Gaussians 
which are slightly shifted from each other. So obviously 
it will show a doM^/e-peaked structure. The locations of 
the extrema of V{q) are easily found in the limit of small 
(j. Note that we are considering a spatially well-localized 
atomic wave packet (Ax <C A) to which the small (3 limit 
applies. Then the locations of the two bumps in ^^{q) 
are given by q « ±V2(Ag) {\-q\ < Ag). Thus the initial 
momentum spread somehow determines the locations of 
the peaks in the split wave packet in the early stage of 
atom-field interaction. 

FIGs.|^(a) and (b) show the shapes of the wave packet 
sitting on a node right after an atomic jump arranged to 
occur at a few different times in the early stage of atom- 
field interaction: (a) in the momentum space and (b) 
in the coordinate space, at various times of the atomic 
jump, r = 1, 2, • • •, 5, 10 and 20 with the same set of pa- 
rameters as in FIG.^. The Gaussian curve is the initial 
packet. Note that the curves are independent trajecto- 
ries. The figures show that our simple argument reaching 
Eq. (|2^) is quite valid up to r ~ 5. It is seen that in the 
coordinate space, the longer the coherent evolution time 
elapsed before the jump, the narrower both the widths 
of, and the separation between, the split peaks become, 
whereas in the momentum space, the wider both become. 

In fact, one can extract two separate dynamical mecha- 
nisms at work in these pictures. Firstly, the packet split- 
ting in the coordinate space — or in other words, the fact 
that V{^) = at a node — right after an atomic jump has 
the following physical ground. The event of an atomic 
jump implies that the atom was in the excited state be- 
fore the jump. For the atom to be lifted from its initial 
ground state to the excited state, the atom-field interac- 
tion must be nonzero. So an event of the atomic jump 
tells us that the probability for the atom to be at the node 
was indeed zero, and our initial Gaussian wave packet has 
been accordingly modified through the back-action of the 
measurement. Secondly, the reason that the separation 
of the two split peaks in the coordinate space becomes 
smaller and smaller as the advent of the atomic jump 
is delayed is explained in the same context. That is, if 
the atomic jump has not occurred after all, we have to 
conclude that actually the atom has not gone through 



n 




(28) 



the Rabi interaction with the cavity field, and has stayed 
in the ground state. This means that the wave packet 
becomes narrower and narrower around a node as time 
goes by before an atomic jump really occurs. When a 
jump occurs, the packet is split into two through the first 
mechanism, but since the separation of the split packet is 
determined by the width of the packet before the jump, 
it is narrower also. Thus one may add that the theory of 
the quantum trajectory is none other than the process of 
continuous correction of our "lack of knowledge" about 
the state of the system, based upon the information pro- 
vided from the measurements. The measurements here 
are of course made by the environment — i.e., the vacuum 
field which is continuously monitoring the atomic damp- 
ing. 

The irregular shifts of the packets in the momentum 
space are due to the random momentum recoils from the 
atomic spontaneous emissions, whereas no such shifts are 
shown in the coordinate space. Note, however, that the 
latter is not because the atomic mass is infinite, but sim- 
ply because 



regardless of the size of the recoil momentum pi or the 
atomic mass, for the dual wave functions ^(a;) and ^p{p). 
The physical reason of this is that the atomic jump has 
been assimicd to occur inst,antaneously, i.e., no time has 
elapsed before and after the jump. So there is no dis- 
placement change before and after the momentum kick. 
However, the absence of motion in the coordinate space 
between the atomic jumps is due to the infinity of the 
atomic mass. As a matter of course, when the atomic 
mass is finite, the system will evolve differently in the 
coordinate space after an atomic jump for a different size 
of the momentum kick that it has gained. Anyhow, we 
see that there is some "internal-motion," if we may, i.e., 
splitting, spreading and separating, of the wave packet 
even for an infinitely massive atom in an optical cavity 
although there is no center-of-mass motion. This inter- 
esting feature seems to deserve some further discussion. 
The following subsection is devoted to this purpose. 

A. The "dynamics" of an infinitely heavy atom. 

Notice that this wave packet splitting in the coordinate 
space by the atomic decays have occurred even for an in- 
finitely heavy atom. Normally, in the usual concept of the 
wave packet "dispersion" introduced in the textbooks, 
the mass of the particle matters. That is, if an infinitely 
heavy particle is in a free space, its wave packet will show 
no dynamics no centcr-of-mass motion or any distortion 
of its shape at all, either in the momentum or in the co- 
ordinate space — as a matter of principle. For the present 
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cavity-QED system with an infinitely heavy atom, how- 
ever, some interesting features appear. The atomic wave 
packet of an infinite mass does show dynamics in the mo- 
mentum space, i.e., sphtting, spreading and separating, 
which is of course due to the momentum exchanges with 
the field. But in the coordinate space, nothing happens 
provided that the system is free of damping, no matter 
what happens to the wave packet in the momentum space, 
as shown in FIGs.|^(a) and (b). In the figure, the inset is 
a lateral view of the 3D-plot on the bottom from the time 
axis which clearly shows the time-invariance of P(^), the 
atomic wave packet in the coordinate space. In the me- 
chanical point of view, this seems rather understandable 
considering that the atomic mass is assumed to be infi- 
nite, but is still quite interesting in the sense that V{^) is 
essentially (but not exactly in fact) a Fourier-type trans- 
form of 'P{q). So, if P(g) changes so radically as shown 
in FIG. 0(a), one would naively expect some change in 
V{(,) also, whatever it may be. 

If we write the wave packet in the same fashion as 
Eq. dl^), V{^) at time t is given by 



(30) 



where a labels the atomic internal states and Vq a nor- 
malizing constant. The coefficients Ca{n,q) continuously 
evolve in time with the Hamiltonian, Eq. dla), but only in 
such a fashion that V{^) remains invariant in time. This 
is in fact easily proved by looking at the time derivative 
ofP(0, i.e.. 



_d 

dt 



ViO^VoJ^J dqdq'e 



iaix(q-q')e. 



Ca{n, q)C:in, q') + Ca{n, q)C:{n, q') 



Now from the dynamical equations, Eqs. (17a) and ( 17b ) 
without the damping term, we have 



Ce{n,q)C:{n,q') = - ^ [C<,(n, g - l)C:(n, g') + C<,(n, g + l)C:(n, g') 



Cg{n,q)C;{n,q') = +^ [C<,(n, g)C:(n, q' - 1) + C<,(n, g)C:(n, q' + 1)] 



etc., and simply because 



j dgdg'e*2-(9-9')« 



-Cg(n, q - l)C:(n, q') + Cg{n, q)C:{n, g' + 1)] = , 



etc., we find 



(34) 



Nevertheless, even more interesting things happen 
when the atom is subject to damping. Then the wave 
packet of an infinitely massive atom does show dynam- 
ical evolutions in the coordinate space as well via the 
coherent decay during the continuous evolution and by 
the atomic jumps. FIGs.p|(a) and (b) show a typical 
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time evolution of a wave packet which has gone through 
a single atomic jump in the momentum space (a) and in 
the coordinate space (b), where 7 — 0.1, Ag — 10, and 
M = 00 for an atom initially in the ground state and 
the field in a one-quantum coherent state in a lossless 
cavity. In FIG.||(b), we clearly see the actual evolution 
of the wave packet in the coordinate space during the 
coherent process, not only the drastic change after the 
atomic jump: The wave packet gets indeed continuously 
squeezed around its peak. The question is how such a mo- 
tion of the wave packet of an infinitely heavy atom can 
occur in the coordinate space. The reason that damp- 
ing causes this type of wave packet evolution even for 



an in finitely massive atom can be seen in Eqs. (17a) and 



( 17t ) . There the atomic damping makes the two equa- 



tions asymmetric in t he sense that it adds a damping 



term only in Eq. ( 17a ) , and thereby we no longer have 
the identities such as Eq. (|33|). 

Since everything happened for an infinitively heavy 
atom, the dynamics following an atomic jump is not a 
mechanical process. That is, the mean position of a 
peak in the coordinate space cannot be related to a time 
integral of that in the momentum space in the fashion 
X = J{p/M)dt, as we did in Sect||, for an obvious rea- 
son. The alert readers will find that the correct relation 
between the mean positions x and p in both spaces right 
after an atomic jump is indeed obtained from the uncer- 
tainty relation, i.e., x ^ a/p where a is a constant on 
the order of h, though it is not exactly h,/2 because the 
wave packet no longer retains the minimum uncertainty, 
as FIGs. ||(a) and (b) demonstrate. Thus, it should be en- 
tirely attributed to the intrinsic wave nature of a matter 
particle as well as to the issue of the position and mo- 
mentum uncertainty which forms the very heart of the 
quantum mechanics p4[]. 



B. Over-damped case 

When the atomic damping is strong such that 7 > 1 
(out of the strong coupling regime), the wave packet es- 
sentially follows the coherent decay only, without quan- 
tum jumps. This can be quickly appreciated from the dy- 



nami cal e quations in the quantum trajectory, Eqs. (24a) 
and ( |24bD . The form of Eq. ( |24a| ) shows that the atomic 
upper-state probability \Ce{n,q)\'^ has an early time de- 
pendence roughly of the form exp(— 7r). FIG.^ shows 
the temporal behaviors of Ve{q) = J2n |C'e(", q)P and 
Pe — J dqV{q) in a trajectory in which an atomic jump 
has not occurred yet. For a large 7, the exponential- type 
factor dictates so that the atomic upper-state probabil- 
ity quickly decays to zero, having no time to grow up to 
an appreciable value. Thus no atomic jumps are likely, 
and the system follows only the coherent decay. While 
damping tends to put the system back to its initial inter- 
nal energy state in this fashion, the atom-field momen- 
tum exchange process rapidly broadens the wave packet 
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in the momentum space, which ends up narrowing — i.e., 
locaUzing — the wave packet in the coordinate space al- 
though the uncertainty product AxAp slowly grows in 
general. This is none other than the phenomenon of the 
"damping-induced localization" (DIL) . For a moder- 
ate atomic damping, the wave packet follows the coherent 
evolution showing the DIL until an atomic jump which 
brings about a sudden splitting — a gigantic change — in 
the wave packet as shown in FIGs.^(a) and (b). 

C. Successive atomic jumps 

The effect of successive atomic decays are studied 
through a numerical calculation, the results of which are 
displayed in FIGs.|l^ (a) and (b), we set 7 = 0.5, k = 
and Aq ~ 10 and M — 00 with the initial field in a one- 
quantum coherent state \a = 1). To quickly appreciate 
the dynamics, we let the system undergo spontaneous 
decays at r = 5, 10, 15, etc. It is shown that the first 
atomic decay splits the momentum space wave packet 
into two peaks, the second into three, the third into four, 
and so on. However, the coordinate space wave packet 
remains largely as a double-peaked shape ever since the 
first atomic decay as shown in FIG.|l^ (b). 

D. Effect of cavity decay 

The cavity quantum jumps, on the other hand, do not 
do much in the wave packet splitting. To see this let us 
turn off the atomic damping channel and open the cavity 
decay channel. Then Eq. ( |2q ) is replaced by the process 

|§) — > a\^) (35) 

which brings about the changes in the coefhcients such 
that 

CJ (n, g) ^ xAICJ (n + 1 , g) (36a) 
C];{n,q) ^ V^^mCl{n + l,q). (36b) 

It will certainly change the shape of the wave packet 
slightly, but never as much as splitting the packet into 
two pieces. Particularly when the mean cavity photon 
number is high, there will be no essential change, whereas 
the packet splitting due to the atomic decay is regardless 
of the intensity of the field. 

IV. EXPERIMENTAL CONSIDERATION 

Since the atomic damping changes the wave packet in 
such a radical fashion, one may expect that it can play 
a significant role in the actual experiments. It is there- 
fore meaningful to compare the result of a hypothetical 
experiment in which the atoms assumed to be nonra- 
diative with that of a realistic experiment in which the 
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atoms are radiative. The actual experimental setup will 
be most possibly such that a beam of atoms are launched 
to fly through the cavity mode while the cavity is con- 
tinuously pumped by an external field, and one measures 
the position distribution of the atoms emerging from the 
interaction region on a surface at some distance from the 
cavity. Since there will be always some fraction of atoms 
which will pass nonradiatively through the interaction 
region without undergoing any jumps even in the pres- 
ence of nonzero damping strength, one can expect some 
nonzero portion of the undeflected atoms in the resulting 
pattern of the position distribution. However, if one can 
still observe two or more split bumps in the background, 
one might say, it is those atoms which have undergone 
any number of jumps that are responsible for the ob- 
served wave packet splitting, particularly considering the 
fact that the wave packet splitting through the coherent 
process tends to be suppressed (by the DIL aforemen- 
tioned, which will be seen in FIGs.|ll] and [ij shortly) in 
the presence of appreciable damping. 

Since it is difficult to analytically discuss the outcome 
of such an experiment, we have performed Monte Carlo 
simulations based on the quantum trajectory theory. As 
we did in our previous work on the cavity-QED atom 
detection system we assume a single mode cavity 
resonant to the transition line of barium (^■^^Ba) atom, 
A = 553 nm, with go/^n = 42 MHz, which correspond to 
/i ~ 1.1 X 10^^. The transverse cavity mode is assumed to 
have a profile so that (7(2?") = (7oexp[— (2/^-|-z^)/it;^] cosfcx, 
where the cavity axis lies in the a;-direction, and w, the 
mode waist is taken to be 37 /im. The atomic longitu- 
dinal velocity is set Vz — 400m/sec. The driving field 
amplitude £ is such that £ / k — 1 so that the mean in- 
tracavity photon number is just one at steady state, i.e., 
{a^a)ss = 1, before the entry of an atom. We choose the 
damping strengths as large as 7/50 — 2k/5o — 0.5. We 
assume that the initial atomic wave packet, to be deter- 
mined by the geometry of the system, is a Gaussian hav- 
ing the momentum uncertainty Ap = IQfik which again 
corresponds to Aa; < A/100, a well-localized, particle- 
like wave packet. Note that the "initial time" is the time 
when the atoms reach the plane normal to the atomic 
path, assumed to be at lOw away, for instance, from the 
cavity axis. Thus, until the wave packet reaches the re- 
gion of appreciable atom-field coupling strength, it will 
freely evolve. 

FIG.|l^ shows the "far- field" position distributions av- 
eraged over 200 atoms detected on a screen located at 
a distance — approximately AQw away, for instance, from 
the cavity axis. (Note that the far-field position distri- 
butions are exactly the same in shape as the momentum 
distributions with proper axis labeling and scaling, ac- 
cording to the diffraction theory.) The thinner line is 
for the case (a) in which the atoms are assumed to be 
nonradiative, while the thicker is for the case (b) where 
the atomic decay channel is let open with 7/30 = 0.5. It 
turns out that during the entire flight time, each atom 
has gone through about 4.4 jumps on average in case (b). 
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The figure shows that, even for the atoms and the cav- 
ity as strongly damped as the given strengths, the wave 
packet splitting robustly shows up, having hardly been 
washed out. The overall width of V{0 {A£_ - 0.68) in 
case (a) is greater than the value (A^ ~ 0.46) in case (b), 
the height of 7^(0) in case (a) being smaller than that in 
case (b) . The narrower distribution of case (b) can be in- 
terpreted as weaker diffusive motion of the atoms in the 
field during the interaction time. Less diffusion implies 
less frequent exchanges of quanta between the atom and 
field, i.e., weaker interaction. So the mean field intensity 
is smaller in the latter case, which again indicates greater 
portion of the vacuum state in the cavity field on aver- 
age. The stronger contamination of the cavity field by 
the vacuum state is due to the introduction of the addi- 
tional damping channel, i.e., the atomic damping, to the 
system. The narrower width in case (b) can be under- 
stood in this way which may be another interpretation of 
DIL. We suggest that the effect of the atomic damping 
will show up as these differences in the shape of the po- 
sition distribution of the detected atoms. Most notably, 
the width of V{^) will be generally narrower than that 
predicted for the ideal systems without damping. 

Let us finally address that the resultant wave packet 
splitting can be entirely attributed to the discontinuous 
atomic jumps, in the viewpoint of the quantum trajec- 
tory theory. For this, we consider a fictitious situation 
(c) where we artificially suppressed the atomic jumps, al- 
lowing only the continuous decay. The result is shown in 
FIG.|lJ by the thin single-peaked line which is in com- 
parison with the lower heavy one, the same curve for the 
realistic case (b) shown in FIG.|l^. As seen in the figure, 
the coherent process alone did not cause splitting of the 
wave packet. On the contrary, it has actually prohibited 
the spreading of the wave packet in the coordinate space. 
(The wave packet would have been much wider than both 
curves even in the free space at the same elapsed time.) 
So it appears that what the atomic jumps did is sim- 
ply to curtail the central portion of the distribution in 
(c) and relocate it to both sides, yielding the double- 
peaked shape in (b). In this regard, we conclude that 
the wave packet splitting in the case of an appreciably 
strong atomic damping is essentially due to the discon- 
tinuous process — in the viewpoint of the quantum tra- 
jectory theory — rather than the coherent process which 
the OSGE was originally envisaged to be. Note also that, 
right after an atomic jump, the atomic state is put into 
the ground state which is simply a 50-50 mixture of the 
orthogonal dressed states. Furthermore, in the presence 
of damping, the final state of the system will eventually 
decay into the product of the atomic ground state and 
the initial coherent field state such that 

lim |*(t)) = l^f) (g> l^a) ^\a)(g> f dpC{p)\p,g) (37) 

where a — £ / n, since the atom-field entanglement gets 
lost as the cavity field is restored to the initial state by 
the driving field whereas the atom eventually decays into 
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the ground state. Thus the atomic states on both peaks 
in the final distribution are again the same. Note that 
without the atomic damping — as the ideal theories usu- 
ally assume, the outcome of such an experiment would 
yield the wave packet split into a pair of distinct bumps 
at large, like the thin line in FIG.|ll|, and each bump 
will correspond to the aggregate of atoms in either of the 
orthogonal dressed-states. When the atomic damping is 
present, however, it is difficult to say that the state on 
one peak is, or has been, orthogonal to that on the other 
because the atomic jumps have been constantly putting 
the atoms on both peaks in the same ground state. 

V. SUMMARY 

We have analyzed the optical Stern-Gerlach effect in 
the realm of the cavity-QED. We included various damp- 
ing mechanisms in the system, and discussed particularly 
the effect the atomic damping in the wave packet evolu- 
tion in the framework of the quantum trajectory theory. 
We have shown that the continuous atomic decay process 
tends to localize the atomic wave packet, keeping it from 
spreading so fast as in the case of no atomic damping, 
whereas the discontinuous atomic jumps tend to split the 
packet. Proposing a possible scheme of an experiment, 
we pointed out that the phenomenon of wave packet split- 
ting could be experimentally observed even in a moderate 
coupling regime where the atomic damping rate is com- 
parable to the atom-field coupling constant, producing 
discernible differences from the predictions made by ideal 
theories in which the damping effects are disregarded. 
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FIG. 1. The model. 

FIG. 2. The evolution of an atomic wave packet placed on a 

node of the cavity field, (a) In the momentum space, and (b) 
in the coordinate space, q = p/hk, ^ = x/X and r = got. For 
a typical atom with fj, = hk^ /2Mgo ~ 1.7 x 10"'*. The atom 
is initially in the ground state having a momentum spread 
Ap = lOhk while the field is prepared in a single quantum 
Fock state. The inset is a lateral view of (b) into the time 
axis along with the cos(fca:) curve. 
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FIG. 3. Spatial modulation of the dressed-energies {h = 1) 
and motion of the dressed-states starting from a node (a) 
and from an antinode (b). "A" and "TV" respectively denote 
antinode and node. 

FIG. 4. Spatial modulation of the dressed-energies {h — 1) 
in the presence of nonzero atom- field detuning, ui' = luq + 5/2 
and g' = ^^fTTW- 

FIG. 5. Wave packet evolution in the coordinate space in 
the presence of nonzero atom-field detuning, 5/go = 0.1, while 
other parameters are the same as in FIG.Q. ^ = x/X. 

FIG. 6. The shapes of the wave packet sitting on a node 
right after an atomic jump arranged to occur at r = 1, ■ • ■, 5, 
10 and 20. q = p/hk, ^ = x/X. 

FIG. 7. The time evolution of the atomic wave packet of 
an infinitely heavy atom free of damping located at a node 
(a) in the momentum space and (b) in the coordinate space. 
The insets are the ground-level views into the time axis. 

FIG. 8. The dynamics of a wave packet in the presence of 
damping in which an atomic jump has occurred, (a) in the mo- 
mentum space and (b) in the coordinate space. 7 — 0.1 and 
Aq — 10, for an infinitely heavy atom initially in the ground 
state placed at a node, with the cavity field in a one-quantum 
coherent state. The insets are the ground-level views along 
the time axis. 

FIG. 9. The evolution of = |Ce(n, Inset 

shows the temporal behavior of Pe = / dqP{q). 

FIG. 10. The shapes of the wave packet sitting on a node 
right after an atomic jump arranged to occur at r = 1, • • ■, 5, 
10 and 20. q = p/hk, ^ = x/X. . 

FIG. 11. The "far-field" position distributions averaged 
over 200 atoms at a location approximately SOw away from 
the cavity axis. The thinner line is for the case of nonradia- 
tive atoms, while the thicker line is for the case of radiative 
atoms with "//go ~ 0.5. 

FIG. 12. The shape of the position distribution that would 
have turned up if all the radiative atoms had gone through 
only the coherent decay (thinner line), compared to the 
thicker line in FIG. 10. Each line averaged over 200 atoms. 
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